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We present the power spectrum of a sinusoidal carrier, frequency 
modulated by a random baseband pulse train in which the signaling-pulse 
duration is finite and the signal pulses may overlap and have different 
shapes. Symbols transmitted during different time slots are assumed to 
be statistically independent and identically distributed. The spectral 
density appears as a Hermitian form suitable for numerical computation 
by a digital computer. Simple conditions in terms of the modulation 
parameters are given under which discrete spectral lines are present in the 
spectrum. Several examples are given to illustrate the method. 

I. INTRODUCTION 

In recent years, digital-frequency and phase-modulation techniques 
have been increasingly important in radio, waveguide, and optical 
communication systems. 

An important parameter in the statistical description of a signal is 
its spectral density, which defines the average power density of the 
signal as a function of frequency. In addition to furnishing an estimate 
of bandwidth requirements, the knowledge of the spectral density is 
also essential in the evaluation of mutual interference between 
channels. 

In this paper, we extend the techniques developed in Ref. 1 for 
digital psk to the case of digital fsk with phase-continuous transitions, 
such as may be obtained at the output of a voltage-controlled oscillator 
driven by a digital baseband wave. We assume that the sinusoidal 
carrier is frequency modulated by a random, baseband pulse train in 
which the signaling pulse duration is finite and the signal pulses may 
overlap and have different shapes. It is generally assumed that the 
symbols transmitted during different time slots are statistically inde- 
pendent and identically distributed. 

We express the spectral density of such an ensemble of continuous- 
phase, constant-envelope, digital, FM waves as a compact Hermitian 
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form that provides an appropriate division between analysis and 
machine computation. The present work permits simpler numerical 
computation of digital fsk spectra than earlier studies, 2-5 and contains 
a simpler statement of the conditions, in terms of the modulation 
parameters, that determine whether discrete spectral lines are present 
in the spectrum or not. 

Examples give the spectra of binary and quaternary fsk waves with 
overlapping baseband modulation pulses of several shapes. 

II. M-ARY FREQUENCY-MODULATED SIGNALS 

We seek the spectrum of the digital frequency-modulated wave : 

x(t) = cos [2irf c t + 0(0], fc > 0. (1) 

*(*) = J'fMdv. (2) 

/«(0- t h, k (t-kT), s k = 1,2, ..-,M. (3) 

*=-co 

The symbol </> is the phase and f d the frequency deviation of the carrier at 
frequency f e . The signaling alphabet consists of M waveforms hi, 
hi, • • • , h Mt that may have different shapes ; one of these is transmitted 
for each signaling interval of duration T. The different signaling 
waveforms in (3) may overlap, but are statistically independent in 
most of the present work; i.e., s k is statistically independent of s t for 
k * I. 

Define for convenience 

j(«) = e»'»i (4) 

then 

x(t) = Re \e^Mt)\- (5) 

The spectral density of v(t) is 

PAD = r $ v (T)e-fi*"dT, (6) 

J —oo 

where 

*,(t) = *^T) s lim ^ f A *,«, r)dt, (7)5 

.4— .oo 6ft J —A 

*.(*, r) = (v(t + r)v*(t)) = < e >l*«+*>-*«>l) 

= (ex V (j J*' fMd»)y (8) 

5 The symbol denotes average on t throughout. 
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An arbitrary constant of integration is implicit in the notation of 
(2), which defines the phase. This is of little consequence in the re- 
mainder of the present paper, since most of our study is directed to- 
ward the spectrum of the complex wave v(t) ; as seen by the final line 
of (8), the absolute phase is irrelevant in determining P v {f)- The 
absolute phase must be rendered explicit only for the following three 
purposes in the present work : 

(i) Relating P x (f) to P r {f) (the spectra of the real fsk wave z(t) 
of (1) and of the complex fsk wave of (4), respectively). 

(ii) Separating possible line-frequency components from the fsk 
wave. 

(iii) Specializing the present fsk treatment to the prior psk results. 1 

To make the absolute phase explicit, we write 

0(0 = f ' fMdn + 0(0) ; »(0) = «*»>. (9) 

Jo 

The term 0(0) is the phase of the fsk wave at t = 0. We consider three 
representative assumptions for (0) : 

(i) 0(0) deterministic, e.g.: 

0(0) = 0. (10) 

(ii) 0(0) random, uniform, and independent of the modulation Sk'. 

Pr [0 < 0(0) £ + d0] = tt> ^ < 2tt. (11) 

(Hi) 0(0) dependent only on the modulation parameters (or signal- 
ing pulses) that contribute to /d(0+) 5 in (3) : 

0(0)= Z ["* h. k (n)dn. (12)§ 

The signal-pulse duration is assumed finite; consequently, in 
(12), the £ has a finite number of terms, and the lower limit 
on the f becomes finite. Specific examples of (12) appear in 
Section IV below. In other words, the net phase contributed by 
all past signaling pulses that are over by t = 0, i.e., for which 
h. t (-kT+) = 0,5 is normalized to 0. 



§ The +'s indicate that / rf and ft, are to be evaluated an infinitesimal time later 
than and —kT, respectively, if any of the signal pulses have discontinuities at their 
upper limits at a time-slot boundary (see Section IV). If the signal pulses and, hence, 
the instantaneous frequency deviation /,/(/) are continuous, these two +'s may be 
dropped. 
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In the appendix we show that, except for special modulations with 
low carrier frequencies f e that take on special values, the following 
simple relation gives the spectrum of the real wave x(t) of (1) : 

PAf) = \PvU-fc)+\P V (-f-fc). (13)5 

The first term of (13) is the spectrum of the complex baseband wave 
i/(0 shifted to the carrier frequency + / e ; the second term is the spec- 
trum of v*(t) shifted to -/ c . More specifically, (13) is valid if any of 
the following are true : 

(i) Eq. (11) holds, independently of any other considerations. 

(u) Eq. (10) or (12) holds, and P„(/) [and hence P,(/)] has no 

line components, 
(m) Eq. (10) or (12) holds, P v (f) [and hence P x (/)] has (equally 
spaced) line components, and f c is low enough so that the two 
terms of (13) overlap, but f c is such that the line components 
of the two terms of (13) do not coincide. The discrete values 
forbidden to f c under these conditions are given by (146) in 
the appendix. 

(iv) The carrier frequency is so large that the two terms of (13) do 
not significantly overlap, independently of any other 
considerations. 

Consequently, we study only P v (f) throughout the remainder of this 
paper. This suffices for all cases except that of a low carrier frequency 
f c that takes on special discrete values related to the baud rate 1/T 
and the modulation, for very special modulations that result in line- 
spectral components (psk is one such case, but there are others). 

The condition expressed by (12) is assumed in much of what follows. 
This results in no significant loss in generality in treating the spectrum 
of the real wave x{t), as discussed above. It permits economy of nota- 
tion in the general case, a convenient treatment of line components 
present in special cases, and simple specialization of the present fsk 
results to prior psk results. 1 

III. NOTATION AND STATISTICAL MODEL 

We introduce the vector-matrix notation of the prior psk study. 1 
Since that study allowed correlated modulation parameters s k while the 
present work is largely restricted to independent s k , only a portion of 
Section III of Ref. 1 need be summarized here. 



§ Compare (5) of Ref. 1. 
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We write (3) as 






where 



& = a* J « W aj 2 > • • • a£*>], (15) 

,h(Q, = h(0 J - [*i(0 MO • • • MO]. (16) 

For a given A: (i.e., for a given time slot) one of the a*'s is unity and the 
rest are zero : 

a («) = i ; a jp = o, i * s k . (17) 

Thus, a* is a unit basis vector, i.e., a* has one component unity and all 
other components zero. 

The modulation process s k is assumed stationary, as in the prior 
psk study, 1 but here we make the stronger assumption of independence 
in most of what follows. Define the first-order probability 

Wi m Pr \s k = i) (18) 

as the probability that the ith signaling waveform is transmitted in 
the A;th time slot. Wi is independent of k by stationarity. By 
independence, 

Pr [sk = i, si = j} = WiWj, k 9* I. (19) 

Then 



Wi*= Pr a*- [0 ■•■ 1 ••• ] 

12 ••• i-1 i t + 1 ••• M 



(20) 



Normalization of the total probability requires 

£ Wi = 1. (21) 

i-l 

We use the following vector notation for the probabilities : 

^ = w]' = [wi w 2 ■ ■ • wjfj (22) 



5 The following notational conventions are adopted throughout : 

(i) Boldface quantities denote matrices, 
(u) Row and column vectors are distinguished by the additional notation i_ i and 

], respectively. 
(in) Ordinary matrix multiplication is indicated by •, Kronecker matrix products 

by X (see second footnote, page 908, Ref. 1, for properties of Kronecker 

products used throughout the present paper), 
(it;) The transpose of a matrix is indicated by ' . 
(v) The Hermitian transpose of a matrix is indicated by t- 
(vi) Multiple Kronecker products are indicated by TJx ( see footnote, page 911, 

Ref. 1 ) and the Kronecker power is indicated by an integer exponent enclosed 

in square brackets. 
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Further define a vector (here of dimension M) with all elements unity 



as 

i-lJ-Cl 1 •••!]• 

Then (21) may be written 

i-w]= L w J -l]= 1. 

Finally, for convenience later, define the diagonal matrix 



w rf = 



W\ 







10-2 







w M 



Then, from (22) and (23), 

w d -l] = w], l l J -w d = w ; 

Note that 

(a*) = w, 

<a*]-a fc > = w d . 
i_i 

IV. FSK AS A BASEBAND PULSE TRAIN 

We seek an expression of v{t), given by (4), (9), (12), and (14), of 
the form 



(23) 
(24) 

(25) 

(26) 

(27) 
(28) 



v(t) = £ fivtit - kT)2 



(29) 



for signaling pulses of finite duration. Assume the hi(t) of (3) are 
strictly time-limited to an interval KT, as follows : 



h(0] = 0], t ^ L K , t > u K . 



T, K odd. 



2 l ' 



K + l 



T, 



Uk = 



K T 
2 1} 



K even. 
Kodd. 
K even. 



(30)1 



(31) 



L K and U K are respectively the lower and upper limits of the pulses. 
Figure 1 shows portions of f d (t) for four different maximum signal- 
pulse durations; the terms k = -1, 0, 1, 2 of (14) are shown, and for 



$&= 0J is a vector of appropriate dimension (here M) with all elements zero. 
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Fig. 1 — Frequency modulation for different signal pulse durations. Index k is 
shown near peak of each pulse. Also, for simplicity, same signal pulse is shown for 
each k. T = time slot duration or signaling period. KT = maximum signal-pulse 
duration. Note different pulse center location for odd and even K. 



convenience a* has been taken the same for each of these time slots. 
The pulses are positioned along the time slots such that the limits of 
each signal pulse lie on the boundary between adjacent time slots (i.e., 
t — integer • T) ; this results in different definitions for Lr, Uk for K 
even and odd. Since symmetric pulses have been chosen for illustration 
in Fig. 1, their maxima are centered in the time slots for K odd, and 
lie on the time-slot boundaries for K even. Discontinuities are permitted 
at the pulse edges (and elsewhere), but are not present in the example 
of Fig. 1 (and would not normally be present in a mathematical model 
of a physical system) ; discontinuities at the time-slot boundaries are 
restricted by the inequalities at the right of (30). Examine the (0, T~] 
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time slot in Fig. 1 as typical ; then the number of pulses contributing 
to fd(t) at every instant equals K} 

It remains for us to express the pulse shapes r{t) and coefficients c* 
of (29) in terms of the signal pulses h(t) and coefficients a fc of (14). We 
give separate treatments for the cases K = 1 and K = 2, and extend 
these results to general K. The treatment is an extension of that for 
the psk case, given in Section IV of Reference 1. 

4.1 Nonoverlapping pulses: K = 1 

The top portion of Fig. 1 shows digital frequency modulation for 
which the signal pulses in different time slots never overlap; in this 
case from (30) to (31), 

h(f)] = 0], t ^ 0, t > T. (32) 

Define 

exp (j J hi(p)dv) exp tj j^ hMdA 

■■■ exp(j J % Mm)*/*)]' 0<t^T. 

q(0- "' ' (33) 

1 ' & t£ 0, t > T. 

Equations (4), (9), (12), (14), and (33) yield 

0(0) = 0; v(0) = 1. (34) 



v(t) = £ S k a k -<i(t - kT)2, (35) 



where 







n a^.q^Dl * < 0- 

i=k+l ' ' 



S k = 1, fc = 0. (36) 

n^-yq(r)l k>0. 



Comparing (35) and (29), the parameters of the latter are given as 
follows for nonoverlapping signal pulses : 

c fc| = S k bk, b fc = &k ; 

r(0] = q(03i 

£* in (37) is, of course, given by (36). 

5 The same conventions were used in Fig. 1, Ref. 1, for the baseband modulation 
pulses in digital psk. 
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4.2 Overlapping pulses: K = 2 

This case is illustrated in the second portion of Fig. 1. In the (0, 2 1 ] 
time slot, the k — 0, 1 pulses contribute. We have from (30) to (31) 

h(0] = 0], t ^ -T, t> T. (38) 

Define 

exp ( j I hMdnj exp h' / h 2 (n)dn) 

• • • exp (j jT A* (/*)<*/*)] . -r < t£ T. 

,q(0, ■ (39) 5 

& i^ -T, t> T. 

Equations (4), (9), (12), (14), and (39) yield 

0(0) = T A.,(M)d/*. (40) 

J-T 

v(t) = t &te-q(< - *r)]Ma*fi-q(* - (/c + l)T)]}, (41) 

t — oo ' ' 

where 5k remains as given in (36), the same as for the prior K = 1 case. 
Proceeding exactly as in (51) of Ref. 1, (41) above yields (29) with the 
following parameters, when no more than two signal pulses overlap : 

c* = S k b,, b^a^Xa^; ^^^ 
r(0] = q(0]Xq«-r)]; 

Sk is given by (36), and X denotes the Kronecker product [see foot- 
note to (14)]. The term b*, like a*, is a unit basis vector, i.e., it has one 
element unity and the remaining M 2 — 1 elements zero. 1 Note from 
(42) and (39) that 

r(0] -0], t £ 0, t > T. (43) 

Binary fsk offers a simple example of these results, governed by the 
same relations between a* and b* as for binary psk given in (57) of 
Ref. 1. 

4.3 Overlapping pulses: general K 

The general case follows by straightforward extension of the above ; 
see (58) to (62) of Ref. 1. Figure 1 illustrates the frequency modulation 
for K = 3, 4. The modulation pulse restrictions are given by (30) and 



§ Comparing (39) with (33), note that the definition of q(0 is different for different 
A'; q(i) ,* over the same interval in which h(t) may be nonzero. 

SIGNAL POWER SPECTRUM 1103 



(31). Define 

exp Ij J hMdfi) exp (j j h 2 (n)dfij 

• • • exp (j j h M (fi)dnj > L K < t ^ Uk- 

q(0 ■ " (44) 

1 ' Q, t£L K , t> U K . 

0(0) of (12) is 

(iC-3)/2 r-kT 

T / hsMdfi, K odd, JC> 1. 

*--(K-l)/2 JL K 

0(0) = (45) 

(K-2)/2 r-kT 

£ / h.Mdp, K even, # > 0. 

<fc--(tf-2)/2 JL K 

The parameters of (29) are 

yt-u/a 
IIx ,**+y X odd. 

i— (Jt-D/2 1 

& - Sife; & = (46) 

"^ K/2 

IIx at+v K even. 

t — (K-2)/2' 
(K-ll/2 

nx q(*-»D], ^odd. 
< — (jc-i)/a 

r(fl] - (47) 

IIx q(<-*r)], Keven. 

i <*-2)/2 

IIx denotes a multiple Kronecker product [see footnote to (14)]. 5* 
is given by 

II a,-- ( g +1)/2; q*(EM], A; < 0; 
i-t+i 

& = 1, A; = 0; K odd. (48) 

n ,a.--(g + i)/v q(Cfo)l fc>0. 

t— i 

II a l _ W2) ;q*(t//c)], fc < 0; 
t-t+i 



i : : i 



S k = 1, k = 0; Keven. (49) 

n >-<g/n -q(frg)]» * > 0. 
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Note that 

r(Q] = 0], t ^ 0, < > 2 1 . (50) 

4.4 Discussion 

It is instructive to obtain the prior psk results 1 by specializing the 
present fsk results. Define the phase shift produced by each signaling 
pulse as 

gi {t) m [' hMdfx, i=l,2,---,M (51) 

Jl k 

or in vector notation 

((0]- f'hW]*, (52) 

and substitute into (44)5 ( or into (33) or (39) for special cases K = 1, 
2). Then the present results of (29), (31), and Section 4.3 are similar to 
the former psk results of (43) and Section 4.3 of Ref. 1, except for the 
factor Sk] in particular, the equations for q(i), b fc , and i(t) have an 
identical form. The additional factor Sk present in the fsk case ac- 
counts for the total phase shift introduced by each of the signaling 
pulses. 

To specialize the present fsk results to the psk case, we require the 
total area of each of the present modulation pulses to be zero. Thus, 
(51) and (52) become 

g t (U K ) m f * hi(t)dt = 0, i = 1, 2, • ••, M (53)55 

Ji K 

or 

g(U K )l= / h(0]#-0]. (54)55 

Substituting in (44), 

q(l7*)] = 1]. (55) 

Therefore, (17) yields 

&-q(£/*)] = 1. (56) 

Substituting in (48) to (49), 

S k = 1, all fc, K, (57) 

completing the specialization of the present fsk results to the psk 
case. 1 

In the general fsk case, the results of Sections 4.1 to 4.3 above reduce 
the fsk problem to determining the spectrum of (29). This is ac- 



I Only the range L* < t ^ U k is relevant in (51), (52), since q(i) ■ outside this 
range. 

S5 Equations (30) to (31) and (53), (54) render (51) and (52) zero for t ^ L K , 
t ^ U K , satisfying (58) of Ref. 1. 
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complished by deleting (149), (153), and (163) to (171) of Ref. 1, 
Appendix B. 5 Then the spectral density of (29) above is given by setting 
b -> c in (150) to (152) and (160) to (162) of Ref. 1 as follows: 

P,{f) = ^R(/).-P c (/D-R»(/)]. (58) 

Pc(/) = f e-*""* c (n) (59)» 

n=— oo 

*c(n) = <c fc+ «]-c5 (60) 

R(/)] = R(/)', = f" e-**"t(tndt. (61) 

1 ' J — 00 

In these relations R(/) is the Fourier transform of r(0 of (37), (42), or 
(47), depending on K (i.e., the amount of pulse overlap). $e(n) is de- 
termined from (36) and (37), (36) and (42), or (46) and (48) or (49) 
in Section VI for K — 1, 2, and general K, respectively. 

V. FSK WITH LINE COMPONENTS 

From (21), (22), and (44), 

|w-q(t/*)]| ^ 1. (62) «» 

We show in the present section that equality in (62) corresponds to 
the presence of line components in the fsk spectrum. Conversely, if 
the inequality of (62) is satisfied, we see in Section VI that the fsk 
spectrum contains no line components. 

Assume throughout the remainder of the present section that 

l l w-q(^)]l = 1. (63) 

This yields 



/, 



Vk 



hi(t)dt = 2irfi + integer -2ir, i = 1, 2, • • -, M, (64) 

where 2vfi is defined as the common area (mod 2w) of all of the fsk 
signaling pulses in the line case. For definiteness we take 

-* < fi £ \. (65) 

Equivalently, 

q(ttr)] = a** 11 ~\ < /i £ !• (66) 



§ These deleted portions were relevant to the study of line spectral components of 
psk. 1 The line spectral components of fsk are treated separately in Section V. 
55 See footnote, page 905, Ref. 1. 
555 This follows because |g.(0l = 1- 
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Now substitute (66) into (48) and (49) ; we have 

S k = e*"*/!, all k, K. (67) 

Thus, when (62) is satisfied 

v(t) = £ e#*"* b^fit - kT)]. (68) 

*=-co 

The relation (68) is the same as the psk result of (43) of Ref. 1 
except for the factor e ik2r/l ; as noted following (52), b* and i(t) have 
the same form in fsk and psk. Consequently, when (63) is satisfied, the 
fsk spectrum may be obtained by simple transformation of the psk 
results of Ref. 1. One way this may be done is to rewrite (68) as 

v(t)(r#*V llT)t = L b k -(e-w< ITHt - kT) i(t - kT)J). (69) 

Comparing with (43) of Ref. 1, all the psk results of Ref. 1 apply 
directly to the present case by making the substitutions 

v{i)->v{t)er**«*i™ 

r(0]->e-> 2 '^" r >'r(0]. K ' 

Therefore, when (68) is satisfied, and consequently (68) holds, the 
fsk wave has line spectra. We separate the line and continuous com- 
ponents as 

»(0 = w,(0 + e (fl. (71) 

The line component is given from (70), and (114) 5 of Ref. 1, as 

i „=_«, \ i /j ^ 72 ^ §§ 

P-,C/) = ji |w^-R(/)]| 2 n ± m ^ (f - Vl ^t 1 )' 

The spectral density of the continuous component is found from (70), 
and (69),*« (96), or (116) of Ref. 1 as follows: 

PuU) = |,R(/).-< >-R*(/)]- (73) 



§ Or from (66) and (95) of Ref. 1 for the special cases k = 1, 2. 

§§ See footnote to (14) ; an exponent enclosed in square brackets denotes the 
Kronecker power, i.e., the Kronecker product of a matrix (or vector) with itself the 
indicated number of times. 

§§§ This is equation (69) of Ref. 1 specialized to independent signal pulses, i.e., 
with (32) used for all n^0, yielding (72) or (73), of Ref. 1. If this restriction is not 
imposed, and (70) above is used in the general form of (69) of Ref. 1, we obtain the 
spectrum of digital fsk with line spectra, with nonoverlapping signal pulses having 
arbitrary correlation (rather than being independent, as in the remainder of the 
present paper). This is the only correlated case that can be readily treated by the 
present methods. 
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{ } = Wrf-w]-^, K = 1. (74) 

+ (w_,X w d X w] - wP> • 1 w, m )e-'" 2,r(/3W ' ) 
+ (w] X w d X ffi - w]i»i v£'"i)« +i,,rC/r " ,l) i -K = 2. (75)5 

, _ expression in { } in (116) of , R (m 

1 } " Ref. 1 with fT->fT- U, generai K ' {7b) 

The condition of (63), that has here been shown sufficient, is shown 
in Section VI also to be necessary for line components to be present 
in fsk spectra. This condition has a simple physical interpretation. 
Every signal pulse must introduce the same total phase change in the 
modulated carrier in order to have line components. This phase change 
has been denoted as 2ir/*; the line components appear at frequencies 

±(/c + £-±-^)> »- ••, -1,0,1, .... (77) 

When fi = 0, (73) to (76) show that the fsk spectra are identical to 
the prior psk results. 1 However the wave in this case is not necessarily 
a psk wave. The stronger condition of (53) or (54) is required to have 
a psk wave ; this condition demands the net phase shift introduced by 
every signaling pulse to be zero. A wave with fi = 0, but one or more 
signal pulses with net phase change equal, for example, to ±2v, will 
have a spectrum given by the psk formula, but will not be a psk wave. 

VI. FSK WITH NO LINE COMPONENTS 

Assume throughout this section that the inequality of (62) is 
satisfied : 

fcq(^)]| < 1. (78) 

We demonstrate that under these conditions the fsk spectrum contains 
no lines. The properties of Kronecker products given in Ref. 1 in the 
second footnote, p. 908, and in the first footnote, p. 915, are used 
throughout without further comment. 

6.7 FSK spectrum with no line component*: nonoverlapplng pulses, K = 1 
From (37) and (60), 

«.(n) = <(&+.<SS (■*+»] ■,«*)>■ (79) 



From (36), 



S k+n St = E &4-*(D3i «• > 0. (80) 



§ See footnote to eq. (72). 
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Let us first consider <$ c (0). Since S k S* k = \S k \ 2 = 1, for all k, we 
have from (79) and (28) 

*e(0) = w d . (81) 



Next, (79) and (80) yield 

n-l 



* e (n) = 



n i ^ ± i-q(T)] {q(7 1 )-a i .]}{a t+n ]-a,}, n £ 1, (82) 



where we have split off the first factor of (80). For n = 1, the first 
{ } , containing the H, is dropped. The first { } is independent of the 
remainder of the expression by (19) ; using the fact that ,q(T), -ajfc] is a 
complex number or equivalently a 1 X 1 complex matrix, and using 
(27) and (28), 

*.(n) = ).w; q(D] } "-'(a.+n] ;q(T\- a*] -,a fc ) 

= [w;q(r)])"-Mw]-,q(r),-w d }, nfcl, (83) 

the last step following from the independence of a A+ „ and a fc . 
Finally, taking the Hermitian transpose 5 of (79), 

*t(n) = (QS^aOM-.UH-Q). (84) 

Alternatively, setting n — > — ?i in (79), 

*.(-n) = ((SyS k - n )(*k-n]-&)). (85) 

Since (79) has been shown independent of k, we may set k — > /c + n 
in (85), to yield upon comparison with (84) 

*.(-n) = *|( n ). (86) 

From (58) to (61), (81), (83), and (86), 

P e (/T) = A + At, (87) 

A = i Wd + e-^ r w]-q(r)-w d £ {e-'" 2 " r l w J -q(T):|}»- 1 . (88) 

1 ' n=l 

Because of (78) the geometric series in (88) converges, and 

1 e -/2*/r w ]. q(r) .w d , , 

A = 2 Wd + l-6-^w-q(T)]- (89) 

Finally, the spectrum of the complex fsk wave of (4) is 

pm) = ^ Rca- ( A + At ) - R *c/)i ( 9 °) 

with A given by (89) for nonoverlapping pulses, K = 1. 



5 See footnote to (14). 
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It is clear that the restriction of (78) renders A and, hence, P v (f) 
finite for all /; consequently, there can be no line components in the 
fsk spectrum. If we take the limit as [w/qCT)]! — > 1, and substitute 
(66) into (88), (90) yields directly the appropriate results for the line 
case, i.e., (72) and (73) to (74) of Section V. 

6.2 FSK spectrum with no line components: overlapping pulses, K = 2 

*.(n) = ((S h+n St)(b k -]-b^) (91) 

with 

fe^&Xfe. (92) 

We have noted in Section 4.4. that the present b fc for fsk has identical 
form to the prior b k of Ref. 1 for psk. In Section 4.2, we saw that S k is 
identical for K = 1 and K = 2 ; therefore, (80) applies to the present 
case as well. Comparing the psk analysis of Section VII of Ref. 1, we 
evaluate (91) above by inserting (80) above inside the ( ) in (87) of 
Ref. 1. 

For n = 0, (87) and (90) of Ref. 1 yield 

*e(0) = wJF. (93) 

For n = 1, (87) and (91) of Ref. 1 yield 

*c(l) = <{.q(n-a*]} X&X {a* +1 ] X^J X a fc+2 ]> 

= <{aiTJ,-a»]- 1 a s j X [a^ij-jat+i,} X a* +2 ]>. (94) 

Since a*, a* + i, and a fc+2 are independent, (27) and (28) yield 

*e(l) = {,q(n-Wd} X w d X w]. (95) 

Finally, substituting (80) inside the ( ) of the third line of (87) of 
Ref. 1, 

*.(n) = /rnV + ,vq(r)]| U(D.-«*]} U(D.-Uf J) <«*f ]•■•*! 

X {a fc+ n +f ]-,a fc+ i,}V n ^ 2, (96) 

where we have split off the first two factors of (80). For n = 2, the 
first { J , containing the IL is dropped. Regarding the factors ^(T 1 ), - a fc ] 
and q(r) -a <;+ i] as complex numbers or alternatively as 1 X 1 complex 
matrices, and using the independence of the different a*, (96) yields 

*.(n) = {ffi-q(T)]}»- 2 <{a, + „]- i q(r2-a fc ]- & } 

X ja t+n+1 ]-,q(r),-a t+i ]-,a fc+ ,J) 
= (w-q(7 1 )]}"- 2 {w]- ,q(r) i -w d } ' 2 ', n ^ 2. (97) 
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The remainder of the analysis proceeds as in Section 6.1. The fsk 
spectrum for overlapping pulses with K = 2 is given by 

P.(/) - ^ .R(/). - (A + At)-R*(/)] ( (98) 



where 

[2] 

(99) 



A = \ w ?3 + r*^tam.w,l X w d X w] + { ^Z'r*j ffl im-] 



Pv(f) again contains no spectral lines by (78) ; the limit l,w,-q(D1 — » 1 
again yields the appropriate results of Section V for the line case. 

6.3 FSK spectrum with no line components: overlapping pulses, general K 

Similar analysis yields the generalization to any overlapping signal 
pulses. From (46) and (60), 

* c (n) = «S k+ jSl)(b k l-&)) (100) 

with b fc given by (46) and S k by (48) to (49). Thus, 

II a*-((jg+i)/»)+,-.- q ( — 2~ TJ , k odd, 

S^ n St = 3 ' n > 0. 

n fegffl+i, - q ( f r ) ] ' fc even ' (101) 

Since the present b* have identical form to b* of Ref. 1, (100) is 
evaluated by substituting (101) inside the ( ) of (106) of Ref. 1, and 
making corresponding changes in the remainder of Section VII, Ref. 1. 
The factor 

Ila^qCtf*)] ( 102 ) 5 

is inserted inside the ( ) of (107) of Ref. 1. In the following we recall 
that the factors of (102), aq] = fl/a], may be regarded alternatively 
as complex numbers or as 1 X 1 complex matrices. 

In (109) of Ref. 1, the first factor of the second line is modified as 

<ai> -> ,q(C/*), • <a<] ■ a*) = (j(Uk\ ■ w d ; (103) 

consequently, 

*e(l) = { ,q(^), -w d } X wS K - 1] X w]. (104) 

In (110), Ref. 1, the first two factors of the fourth line are modified 



fi U K is given by eq. (31). 
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as in (103), yielding 

4 e (2) = { ,q(tM, -w d }' 2 ' X w5 x - 2 ^ X wj»1. (105) 

By induction, 

*e(n) = { ,q(t/g), -w d } w x wS*"" 3 X w]«, n ^ K. (106) 

Next, for n ^ K, inserting (102) inside the ( ) of the second line 
of (107) of Ref. 1, associating the factors of (102) with corresponding 
factors of the second II x of (107), Ref. 1, with n - K factors of (102) 
left over, and then noting that all factors have different indices and, 
hence, are independent, 

•c(n) = w]'" • { ,q(C/g), -w d } « l .qd/jr). •*]}-* 

= {■w,-q(g7jr)]}- ic {w]- ,q(yjr), -w 4 »}w > n ^ X. (107) 

Finally, from (108) of Ref. 1 and (100) above, 

*e(0) = wS K] . (108) 

The fsk spectrum for overlapping pulses with general K is given by 

P v (f) = ^ ,R(/), - (A + A+)-R*(/)], (109) 

where 

A = i wS x] + *£ «-**"* <L(Uk)M * 2 X w£*-" 3 X w]W 

+ 1 _ e -^/r Wj . q ([/ x )] ' ^ iU J 

The condition (78) again guarantees no spectral lines; as equality is 
approached in (78) the present results approach those of Section V for 
the line case. The condition (63) is, therefore, necessary and sufficient 
for P v (f) to have line components. 

VII. ILLUSTRATIVE EXAMPLES 

The computation of the digital fm spectral density from the above 
methods is straightforward. For a given set of baseband signaling pulse 
shapes and symbol probability distribution, we determine K, the over- 
lap parameter, and the probability row vector jr, We then evaluate 
q(U K ) and ^-q (£/*)]. 

If \w.-a(Uir)\ < 1, we know that there are no line components in 
P v (f). The continuous part of the spectrum is evaluated from the ap- 
propriate Hermitian form given in Section VI. 

If [w.-q(Uk)1\ = 1, we know that there are line components in the 
spectrum, and also that 

q(t/x)] = e*"M], |/i| £i (HI) 

We determine fi, and then P Vl (f) and P Vc (f) from the methods given 
in Section V. 
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In the following examples, the digital computer is programmed to 
work directly with the Hermitian forms (both ordinary and Kronecker 
matrix multiplications are performed by the computer). In this way, 
complicated cases involving multilevel signal pulses overlapping 
several time slots may be simply treated. 

The case of rectangular-pulse fsk modulation is treated in Ref. 3, 
and consequently will not be considered here. 5 

We make the following assumptions for convenience; none are 
essential. 

(i) The number of frequency levels is a power of 2, 

M = 2 N , N an integer. (112) 

(it) The M baseband signaling pulses have a common shape ; 

g k (t) = C k g(t). (113) 

(in) All signal pulses are equally likely ; 

w] = -^1], w d = jjI M , (114) 

where Iat is the identity matrix of order M. 

7.1 Raised-cosine nonoverlapping signal pulses: K = 1 

If the pulses have a common raised-cosine shape, 



A! 

A 2 



irfd 



1 - cos 2 ~y 1 , < t g T, f d > 0, 



g(0] = 

0], otherwise, (115) 

where Ai, A 2 , • • • , Am are the peak-frequency-deviation parameters of 
the fsk signals. We assume that Ai = 1, A 2 = — 1, A 3 = 3, A 4 = —3, 
-.., A„= -(M - 1). 
Since K = I, 

R k (f) = T £ Jn ( fj ^) e-*"-<"WT 

sin* (fT-f^ + n) 

— r — t^t — \ ' (116) 

where J »(x) is the Bessel function of the first kind and of order n. 

I This is the only discrete-frequency-modulation spectrum given in Ref. 3. 
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7.1.1 Raised-cosine signaling with no line spectrum: K = 1 

Since 

q(£/j) = p&irfdT QJbirfdT . . . QJbMTfdT^ 



(117) 



line spectra are absent if and only if f d T is not an integer. In this case, 



p.U) = Puif) = *£(/),• (A + At).R(/)!, 



T 



(118) 




Fig. 2— Spectral density of binary fsk system with raised-cosine signaling and 
pulse duration T. K = 1. 2f d is the spacing between two adjacent a priori chosen 
frequencies. 



1114 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1975 



? - 







(b) 


1 (0.198) 






_ 


(0.194) 






- 


1 1 ' 


/fdT=1.5 

/ 
/ 
/ 

/ 

/1.75 

/ 

/ 

^1.25 

\S ^ 2.00 


► (0.0481) 






>Jvv/y7 \"\"^^ 


\ 


1 



1.5 2.0 

fT 



2.5 



Fig. 2 (continued). 



where R(/)] is given by (116) and 



1 l e-' w l]-,q(E/i). 

2Af U + M» _^ /T _1_ sin ilf^/,7 1 

6 M smirf d T 



3.0 



(119) 



For M = 2, 4, and 8 and various values of k = /^T, P„ (/) is plotted 
in Figs. 2, 3, and 4. 

7.1.2 Raised-cosine signaling with line spectrum: K — 1 

The fm spectral density P v (/) contains lines if and only if | L w J - q( £/i)] | 
= 1 ; that is, if and only if 

f d T= 1,2,3, ■•■. (120) 
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Fig. 3— Spectral density of quaternary fsk system with raised-cosine signaling 
and pulse duration T.K = 1. 2/ d is the spacing between two adjacent a priori chosen 
frequencies. 
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Fig. 4 — Spectral density of octonary FSK system with raised-cosine signaling and 

fiulse duration T. K = 1.2/j is the spacing between two adjacent a priori chosen 
reqi; 



frequencies. 
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In this case, 



fi — I. SdT — 1, 3, 5, 
Si = 0, f d T = 2, 4, 6, 



(121) 
(122) 



From Section V, 



(*-•-¥) 



+ ^RUl- (B + Bt).R(/)t, (123) 
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Fig. 5 — Spectral density of binary fsk system with raised-cosine signaling and 
pulse duration 2T '. A' = 2.2/,; is the spacing between two adjacent a priori chosen 
frequencies. 
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Fig. 5 (continued). 



where 8 



B + Bt = F I --^ 1 >A- 



(124) 



For M = 2, 4, and 8 and for some integral values of f d T, P v (f) is 
also plotted in Figs. 2, 3, and 4. 



§ Note that B + B* for K = 1 is given in (74). 
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f d T=0.125 



(a) 




(4.064X10" 



.0.375 



(1.192X10 -4 ) 



Fig. 6— Spectral density of quaternary fsk system with raised-cosine signaling 
and pulse duration 2T. K = 2.2f d is the spacing between two adjacent a priori 
chosen frequencies. 
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Fig. 6 (continued). 

7.2 Raised cosine overlapping signal pulses: K = 2 

If a raised-cosine signal pulse just fills up two time slots, 



Ai 

A, 

A.u 



»/-( 



1 + cos ^ ) ■ T < I g 7'. /, > (! 



8(0] - 



0], 



otherwise, (125) 



K = 2, and the spectral density may be calculated from Sections 4.2, 
V, or 6.2 according to whether 

lw,.q(t/ 2 )]| = 1 (126) 



or 



(127) 



Lw/q(( T .>)! < 1. 

If A, = 1, A 2 = -1, A 3 = 3, A 4 = -3, • • •, Ajr-i = (M - 1), 
A.i, = — (M — 1), (116) can be shown to be satisfied if and only if 
2f/r is an integer. Further, 

(128) 



and 



/, = 0, iifdT = 1,2,3, •••, 

/, - h if/-r= if,!, •••• (129) 
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In this case, note that 

&^ T ^ ~ M sin (2»/,D ' ( } 

In this case, the fsk spectrum contains line components, and the con- 
tinuous and line spectra for M = 2 and 4 are shown in Figs. 5 and 6. 

If 2f d T is not an integer, the fsk spectrum does not contain any lines, 
and the continuous spectrum given by (97) to (98) for M = 2 and 4 
is also plotted in Figs. 5 and 6. 

Several observations can be made from Figs. 2 to 6. For both 2^=1 
and K = 2, discrete spectral lines appear as f d T approaches the limiting 
value and the power in the lines is substantial for K = 1 and M = 2. 
Also note that power in the line components with K = 2 is smaller 
than the power in the lines with K = 1. 

For the same value of fdT and K, the principal portion of the spec- 
trum of binary fsk is narrower than that of quaternary fsk. For K = 1, 
quaternary fsk spectrum is narrower than that of octonary fsk spec- 
trum for the same value of fdT. 

Since lines can appear in the spectrum for a set of values of f d T, 
the fsk spectrum is quite different from the psk spectrum even when 
lines are present in the fsk spectrum. 

VIII. SUMMARY AND CONCLUSIONS 

Matrix methods are given to express the spectral density of a carrier, 
frequency-modulated by a random baseband pulse train, in a concise 
and computable form. Arbitrary pulse shapes may be used for Af-ary 
digital signaling, and they may overlap over a finite number of signal 
intervals. 

The spectral density is expressed as a compact Hermitian form 
suitable for numerical computation by a digital computer. The com- 
puter is readily programmed to perform matrix operations directly, 
rather than expanding the Hermitian form and evaluating the indi- 
vidual terms. In this way, quite complicated cases involving multilevel 
signal pulses overlapping several time slots may be treated simply. 

Simple conditions in terms of the modulation parameters are given 
under which discrete spectral lines are present in the spectrum. A 
method is given to evaluate the power in the discrete spectral lines. 
The utility of the method is illustrated by giving several examples. 

The present results are restricted to independent signal pulses, 
except for the very special case of nonoverlapping (K = 1) signal 
pulses with line spectra present [see footnote following (72)]. In Ref. 1 
we saw that for overlapping signal pulses with K = 2 (i.e., signal 
pulses that occupy no more than two time slots), psk with correlated 
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signal pulses required the fourth-order statistics of the digital modula- 
tion process, fsk is more complicated; the correlated K = 2 case 
requires statistics of all orders for the modulation process, except in 
the line spectrum case. However, certain special cases with correlated 
signal pulses have been handled by extension of the present methods. 

APPENDIX 

Spectra of Complex and Real FSK Waves 

From (130), (131) of Ref. 1, Appendix A, we see that (13) of the 
present paper holds true if 

e>*"''*„.(*, r) = 0, (131) 

where 

*...«, r) = (v(t + tMO) = <e>i*«+*>+*«>'). (132) 

From (3) and (9), 

*.,.(*, r) = /•*♦« exp (j h jt m [f^ KT + j^ k J\ h.Mdf^y 

(133) 

If 0(0) is independent of the modulation parameters s k , (133) 
becomes 

*,..«, r) = <•**»>) ^exp (j h £ m [j^ kT + /J," J h -Md»)y 

(134) 

Thus, if (11) holds, the first factor in (134) = 0, and (131) and, hence, 
(13) follow immediately. 

If, instead, (10) holds, and imposing in addition the independence of 
the St, (134) becomes 

*..•(«, r) = k fl m (exp (j [ j* T ~ >T + J^ T ] h.Mdn)y (135) 

For large enough |^| and finite pulse length, the integrals in the ex- 
ponent fall into three classes, depending on A; : 

(t) Both limits on both integrals lie to one side of the signal pulses. 
The integrals equal and the corresponding factor in the in- 
finite product equals 1 and, hence, may be ignored. 
(it) The limits on integrals straddle the pulse, and the limits may 
consequently be replaced by f u L K K , where the maximum pulse 
length is KT and L K , U K are given in (31). 
(Hi) The common lower limit, or one or both upper limits, lie within 
the pulse. 
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Therefore, for large enough \t\ (depending on the pulse length and 
on t) 

$„*(< + T, t) = *,„.(*, r) /exp (j2 J K h n (ji)dv\\ . 

t sufficiently positive. (136) 

*.„.(/ - T,t) = *„»(*, r) /exp f -j2 ^ * h, k (n)d^ . 

J sufficiently negative. (137) 

Condition (12) may be substituted for condition (10) with only minor 
changes in the above discussion, and identical results (136) and (137). 
From (44), 



q«'*)] = 



/ r u K 

U J hi(ji)dp 



)i 

exp (j J h%(ji)dnj 
exp (j I h M (ji)dnj 



(138) 



the integrals in the exponents are the areas of the signal pulses. We 
consider the two cases of Sections V and VI : 

(t) All pulse areas are identical (mod 2tt) ; line components are 

present in the spectrum, 
(u) The contrary ; no line components are present. 



In case (i), from (6) 

q(£M] = e' 2 *"-l], -i</i£ *; 



(139) 



2t/i is the common area (mod 2x) of all the signal pulses and 1] is the 
unit column vector of (23). It therefore follows that the final factors of 
(136) and (137) satisfy the following: 

/exp(±j2 f * h. h (p)dp) \ = e±» T ", line spectrum present. (140) 

/exp(±j'2 I * h 8k (fi)dn)\\ < 1, line spectrum absent. (141) 

With line spectra absent, (141), (136), and (137) show immediately 
that (131) holds for all f e , and, hence, condition (10) or (12) guarantees 
the result of (13) for independent s k without further restriction. 

With line components present, (140), (136), and (137) show that 
*•••(<, r) regarded as a function of t contains a periodic (line) com- 
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ponent. Therefore, we write 

*„„.(/, r) = *„,({, t)i + *„„.(*, r) c , (142) 

where 

*» !.*(', t) c = 0, |£| sufficiently large, (143) 

and 

M(,r)i-«* /,wn £ *>«(T)e* 2 *< l/r >. (144)5 

Only the line component (144) can possibly contribute to the average 
on the left-hand side of (131) : 



e*'' e, *„.(«, t) 

= E *».M lim ^j- T 4 e^lv-rHVi-HO/ri^ 

f r \ i im sin2ir[2/.+ (2/ I + n)/r]il , , 

= J? . ^" (r) .1^1 2,[2/. + (2/, + »)/T]A (MB) 

The limAH.. — > if 2/« + (2/, + n)/T ^ for every integer n, thus 
satisfying (131). Hence, condition (10) or (12) and 

2 (/• + $)*?' n = ' 1 ' 2 ' 3 ' ••' (146) 

guarantees the result of (13) for independent S*. Since, by (72), the line 
spectral components of v(t) occur at (ft + n)/T, (146) is equivalent to 
requiring that the line components of v{l)e> 2w/cl and of v*(t)e-' 2wfct 
never coincide. Note that condition (5) of Ref. 1 is the special case of 
(146) above for fi = 0. 

Finally, (135) to (140) of Ref. 1, Appendix A apply also in the present 
case, and establish that (13) is a good approximation if f e is high enough 
so that there is no significant overlap between the two terms of (13). 
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